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I discuss imaging sources from low relative- velocity correlations in heavy-ion reac- 
tions. When the correlation is dominated by interference, we can obtain the images 
by Fourier transforming the correlation function. In the general case, we may use 
the method of optimized discretization. This method stabilizes the inversion by 
adapting the resolution of the source to the experimental error and to the mea- 
sured velocities. The imaged sources contain information on freeze-out density, 
phase-space density, and resonance decays, among other things. 



1 Introduction 

Imaging techniques are used in many diverse areas such as geophysics, astron- 
omy, medical diagnostics, and pohce work. The goal of imaging varies widely 
from determining the density of the Earth's interior to reading license plates 
from blurred photographs to issue speeding fines. My own interest in the prob- 
lem stems from seeing an image of Betelguese, a red giant ~ 600 ly away that 
has irregular features changing with time. The image..was obtained using in- 
tensity interferometry such as used in nuclear physicsu. After seeing this, the 
natural question was whether images could be obtained for nuclear reactions. 
Needless to say, answers to such questions tend to be negative. 

In a typical imaging problem, the measurements yield a function (in our 
case, the correlation function C) which is related in a linear fashion to the 
function of interest (in our case, the source function S): 

C{q)= I drK{q,r)S{r). (1) 

In other words, given the data for C with errors, the task of imaging is the 
determination of the source function S. Generally, this requires an inversion 
of the kernel K. The more singular the kernel K, the better the chances for 
a successful restoration of S. 

In reactions with many particles in the final state, there is a linear rela- 
tion of the type (|^) between the two-particle cross section dPa/cPpi cPp2 and 
the unnormalized relative distribution of emission points S' for two particles. 

"^Talk given at the Int. Workshop on "Collective Excitations in Fermi and Bose Systems", 
Serra Negra, Brazil, September 14-17, 1998. 
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Figure 1: Separation of the interference and final-state interactions, in terms of the two- 
particle wavefunction, from the amplitude for the reaction. 



Interference and interaction terms between the two particles of interest may 
be separated out from the general amplitude for the reaction and described 
in terms of the two-particle wavefunction (see Fig. |l|). The rest of the 

amplitude squared, integrated in the cross section over unobserved particles, 
yields the unnormalized Wigner function S' for the distribution of emission 
points written here in the two-particle frame: 



J d^rS'pir)\^'cl.^{r)\\ (2) 



The vector r is the relative separation between emission points and the equation 
refers to the case of particles with equal masses. The size of the source S' is of 
the order of the spatial extent of reaction. The possibility of probing structures 
of this size arises when the wave-function modulus squared, |<i>*^~^p, possesses 
pronounced structures, either due to interaction or symmetrization, that vary 
rapidly with the relative momentum, typically at low momenta. The two- 
particle cross section can be normalized to the single-particle cross sections to 
yield the correlation function C: 

CiP^ - P^) = = / Spir) . (3) 

d-*pi d^P2 

The source S is normalized to 1 as, for large relative momenta, C is close to 1 
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and |<i>p in (||) averages to 1: 



jd\Sp{f)^l. (4) 

Depending on how the particles are emitted from a reaction, the source 
may have different features. For a prompt emission, we expect the source to be 
compact and generally isotropic. In the case of prolonged emission, we expect 
the source to be elongated along the pair momentum, as the emitting system 
moves in the two-particle cm. Finally, in the case of secondary decays, we 
expect the source may have an extended tail. 

In the following, I shall discuss restoringJjhe sources in heavy-ion reactions 
and extracting information from the imagesH. 



2 Imaging in the Reactions 

The interesting part of the correlation function is its deviation from 1 so we 
rewrite (0) 



(fir 



- 1 Spif) 



(firK{q,r)Sp{r). 



(5) 



From (^), it is apparent that to make the imaging possible <I>ir'(r) must 
deviate from 1 either on account of symmetrization or interaction within the 
pair. The angle-averaged version of (j^) is 



np{q) = AT: J drr^Ko{q,r)SUr) 



(6) 



where Kq is the angle-averaged kernel. 

Let us first take the case of identical bosons with negligible interaction, 
such as neutral pions or gammas. The two-particle wavefunction is then 



(7) 



The interference term causes 



to deviate from 1 and 



K{q, r) — cos (2<f • r). 
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Figure 2: Gamow-corrected vr tt correlation fugetion for Au + Au reaction at 
10.8 GeV/nuclcon obtained by the E877 collaboration Q. 



In this case, the source is an inverse Fourier cosine-transform oiTZp. Also, the 
angle averaged source can be determined from a Fourier transformation (FT) 
of the angle-averaged C as the averaged kernel is 

K%,r)^^J^. (9) 
2qr 

While neutral pion and gamma correlations functions are difficult to mea- 
sure, charged pion correlations functions are not. The charged pion correlations 
are often corrected approximately for the pion Coulomb interactions allowing 
for the use of FT in the pion source determination. In Figure ^, I show one 
such corrected correlation function for negative pions from the Au + Au re- 
action, at 10.8 GeV/nucleon from the measurements by E877 collaboration at 
AGS a. In Figure ^, I show the relative distribution of emission points for 
negative pions obtained through the FT of the correlation function in Fig. ^ 
The FT has been cut off at Qmax = 200 AleV/ c giving a resolution in the source 
of Ar > 1/(2 qmax) = 2.0 fm. The data spacing gives the highest distances that 
can be studied with FT of r^ax H 1/(2 Ag) = 20 fm. As you see, the relative 
source has a roughly Gaussian shape. 
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3 Perils of Inversion 

For many particle pairs, such as proton pairs, interactions cannot be ignored 
and the straightforward FT cannot be used. Indeed, even in the charged-pion 
case, one might want to avoid the approximate Coulomb correction. In lieu of 
this, we can simply discretize the source and find the source that minimizes 
the x^- This procedure could work for any particle pair. 

With measurements of C at relative momenta {qi \ and assuming the source 
is constant over intervals {Arj}, we can write Eq. (|6|) as 

= Coofe)-l = ^47rArr2Xofe,r,)5(r,) (10) 

j 

^ (11) 

j 

The values Sj can be varied to minimize the x^: 
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Figure 4: Two-proton correlatitw function for the ^'*N + ■^^Al reaction at 75 MeV/nucleon 
from the measurements of Ref. □ for three gates of total momentum imposed on protons 
emitted in the vicinity of ^i^b = 25° . 



Derivatives of the with respect to S give hnear algebraic eqs. for S: 



(13) 



with the solution in a schematic matrix form: 



(14) 



There is an issue in the above: how do we discretize the source? The FT 
used before suggests fixed-size bins, e.g. Ar — 2 fm. However fixed size bins 
may not be ideal for all situations as I will illustrate usinft^Fig. ^. This figure 
shows the pp correlation function from the measurements □ of the ^^N + ^^Al 
reaction at 75 MeV/nucleon, in different intervals of total pair momentum. The 
different regions in relative momentum are associated with different physics of 
the correlation function. For example, the peak around q ^ 20 MeV/c is 
associated with the ^5*0 resonance of the wavefunction with a characteristic 
scale of fm - this gives access to a short range structure of the source. On the 
other hand, the decline in the correlation function at low momenta is associated 
with the Coulomb repulsion that dominates at large proton separation and 
gives access to the source up to (20-30) fm or more, depending on how low 
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momenta are available for C . Should we continue at the resolution of Ar > 2 f m 
up to such distances? No! At some point there would not be enough many 
data points to determine the required number of source values! Somehow, we 
should let the resolution vary, depending on the scale at which we look. 

A further issue is that the errors on the source may explode in certain 
cases. The errors are given by the inverse square of the kernel: 

A^S.^iR-^ K)-^. (15) 

The square of the kernel may be diagonalized: 

(K^ = ^i^fc, i^fc, ^Y^Kut u'^ . (16) 

k ^ ct 

where {u"} are orthonormal and Aq > 0; the number of vectors of equals the 
number of r pts. The errors can be expressed as 

Q 

You can see from the last equation that the errors blow up, and inversion 
problem becomes unstable, if one or more A's approach zero. This must happen 
when K maps a region to zero (remember K = |$p — 1), or when K is too 
smooth and/or too high resolution is demanded. A A close to may be also 
hit by accident. 

The stability issue is illustrated with Figs. ^ and |^. Figure ^ shows correla- 
tion functions from model sources with small errors added on. Figure |^ shows 
the source in 7 fixed-size intervals of Ar ~ 2 fm. This source was restored fol- 
lowing Eq. (p^, from the correlation function indicated in Fig. ^. The errors 
in this case far exceed the original source function. Every second value of the 
restored source is negative. 

Vast literature, extending back nearly 75 years, exists on stability in in- 
version. One of the first researchers to recognize the difficulty, Hadamard, in 
1923 B, argued that the potentially unstable pmblems should not be tackled. 
A major step forward was made by TikhonovEl who has shown that placing 
constraints on the solution can have a dramatic stabilizing effect. In determin- 
ing the source from data, we developed a method of optimized discretisation 
for the source which yields stable results even without any constraints El. 

In our method, we first concentrate on the errors. We use the g-values for 
which the correlation function is determined and the errors of {cr,}, but we 
disregard the values {Ci}. We optimize the binning for the source function to 
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Figure 5: The solid line represents the correlation function from a Gaussian model source 
while the dashed lines represent the correlation functions from a source with an extended 
tail.i-jThe points represent values of C with errors that are typical for the measurements in 
Ref.y (Fig. |). 

minimize expected errors relative to a rough guess on the source 5™°'^: 

i 3 j j \ a " I 

Only afterwards we use {Ci\ to determine source values Sj with the optimized 
binning. This consistently yields small errors and an introduction of constraints 
may additionally reduce those errors. The proton source imaged using the 
optimized binning from the correlation function in Fig. [s] is shown in Fig. |^. 

4 pp Sources 

Upon testing the method, we apply it to the 75 MeV/nucleon + ^^Al data 
by Gong et al. a shown in Fig. ^. In terms of the radial wavefunctions g, the 
angle-averaged pp kernel is 

^o(g,r) = i5](2j + l)(gjf(r)f -1. (19) 

jsW 
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Figure 6: The solid histogram is the relative pp source function S restored from the sim- 
ulated correlation function in Fig. H from the Gaussian model source (open symbols there). 
The dashed line is the original source function that we used to generate the correlation 
function. We employed fixed-size intervals of Ar = 2 fm and we imposed no constraints 
on S. 



We calcuJate the wavefunctions by solving radial Schrodinger equations with 
REID93Q and Coulomb potentials. The sources restored in the three total mo- 
mentum intervals are shown in Fie^H, together with sources obtained directly 
from a Boltzmann equation modelB (BEM) for heavy-ion reactions. 

The sources become more focussed around r = as total momentum 
increases. Now, the value of the source as r — > gives information on the 
average density at freeze-out, on space-averaged phase-space density, and on 
the entropy per nuclcon. The freeze-out density may be estimated from 

p/„e^e=iiVpartX5(r-^0), (20) 

where A'part is participant multiplicity. Using the intermediate momentum 
range, we find 

P/ree.e « ( 17) (.0015fm-=^) = .16po • (21) 

The space-averaged phase-space density may be estimated from 
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Figure 7: Relative pp source function S restored^solid histogram) through the optimized 
discretization from the correlation function in Fig. H (open symbols there), together with the 
original source function (dashed line). 



Using the intermediate momentum range we get (/) « .23 for this reaction. 

The transport model reproduces the low-r features of the sources, including 
the increased focusing as the total momentum increases. The average freeze- 
out density obtained directly within the model is p freeze — •14po- Despite 
the agreement at low r between the data and the model, we see important 
discrepancies at large r. I discuss these next. 

An important quantity characterizing images is the portion of the source 
below a certain distance (e.g. the maximum r imaged): 

A(w) - / dV5(r). (23) 

J r<r,nax 

If fmax =^ oo, then A approaches unity. The value of A < 1 signals that some of 
the strength of S lies outside of the imaged region. The imaged region is limited 
in practice by the available information on details of C at very-low q. We can 
expect pronounced effects for secondary decays or for long source lifetimes. If 
some particles stem from decays of long-lived resonances, they may be emitted 
far from any other particles and contribute to 5 at r > r^nax- 

Table gives the integrals of the imaged sources together with the integrals 
of the sources from BEM over the same spatial region. Significant strength 
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Figure 8: Relative source for protons emitted from the + ^^Al reaction at 

75 MeV/nucleon, in the vicinity of fliab = 25°, within three intervals of total momentum 
of 270-390 MeV/c (left panel), 450-780 MeV/c (center panel), and 840-1230 MeV/c (rig 
panel). Solid and dotted lines indicate, respectively, the source values extracted from datal 
and obtained within the Boltzmann-equation calculation. 



P-Range 






771 ax 


[MeV/c] 


restored 


BEM 


[fm] 


270-390 


0.69±0.15 


0.98 


20.0 


450-780 


0.574±0.053 


0.91 


18.8 


840-1230 


0.87±0.14 


0.88 


20.8 



Table 1: Integrals of sources from data and BEM in the three intervals of total momentum. 



is missing from the imaged sources in the low and intermediate momentum 
intervals. BEM agrees with data in the highest momentum interval but not 
in the two lower- momentum intervals. In BEM there is no intermediate mass 
fragment (IMF) production. The IMFs might be produced in excited states 
and, by decaying, contribute protons with low momenta spread out over large 
spatial distances. Information on this possibility can be obtained by examining 
the IMF correlation functions. 
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Figure 9: Fteiative source for IMFs emitted from central *'*Kr + '^^'^Au reactions from the 
data of Ref. Ej at 35 (dotted line), 55 (dashed line), and 70 MeV/nucleon (solid line). The 
insert shows the source multiplied by . In both plots, the full image extends out to 90 fm. 



5 IMF Sources 

Because of the large charges > 3), the kernel in the case of IMFs is domi- 
nated by Coulomb repulsion. With many partial waves contributing, the kernel 
approaches the classical limit 0: 

Xo(g,r) =^^(r-r,)(l-re/r)l/2-l, (24) 

where — 2iiZiZ2e^ / is the distance of closest approach. There are no 
IMF correlation data available for the same reaction used to measure the pp 
correlation data, so we use data within the same beam energy range, i^. the 
s^xr -19^ Au at 35, 55, and 70 MeV/nucleon data by Hamilton et a/.0. The 
extracted relative IMF sources are shown in Fig. |[ The source integrals for 
the IMF sources are given in Table ^. Interestingly, we are nearly capable 
of restoring the complete IMF sources. For the relative distances that are 
accessible using the pp correlations (~ 20 fm) we find only (70-80)% of the 
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Beam Energy 


A(90fm) 


A(20fm) 


[MeV/A] 






35 


0.96±0.07 


0.72±0.04 


55 


0.97±0.06 


0.78±0.03 


70 


0.99±0.05 


0.79±0.03 



Table 2: Comparison of the integrals of the midrapidity IMF source function, \{rmax), in 
central *'*Kr + ^^^Au reactions at three bffi«n energies, for different truncation points, rmax- 
The restored sources use the data of Ref . tj. 



IMF sources. This is is comparable to what we see for the lowest-momentum 
pp source but above the intermediate-momentum proton source. We should 
mention that we can not expect complete quantitative agreement, even if the 
data were from the same reaction and pertained to the same particle- velocity 
range. This is due partly to the fact that more protons than final IMFs can 
stem from secondary decays. 

6 TT^ VS. Sources 

We end our discussion of imaging by presenting sources obtained for pions and 
kaons from central Au -I- Au reactions at about 11 GeV/nuclcon. This time 
we use the optimized discretization technique rather than the combination of 
approximate Coulomb corrections and the FT. For both meson pairs the kernel 
Kq is given by a sum over partial waves: 

^ElM^-i- (^^) 

where gi{r)s stem from solving the radial Klein-Gordon equation with strong 
and Coulomb interactions. In practice the strong interactions had barely any 
effect on the kernels and the extracted sources. ^ 

The data come from the reactions at 10.8 and 11.4 GcV/nuclcond. The 
respective tt^ and K'^ sources are displayed Fig. ^ The kaon source is far 
more compact than the pion source and there are several effects that contribute 
to this difference. First, kaons have lower scattering cross sections than pions, 
making it easier for kaons to leave the system early. Second, fewer kaons than 
pions descend from long-lived resonances. Next, due to their higher mass, the 
average kaon has a lower speed than the average pion, making the kaons less 
sensitive to lifetime effects. Finally, the kaons are more sensitive to collective 
motion than pions, enhancing the kaons' space-momentum correlations. Differ- 
ences, qualitatively similar to those seen in Fig. O, in the spatial distributions 
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Figure 10: Relative sources of 7r~|-(«:ircles) and of (triangles) extracted frowi central 
Au + Au data at 11.4 GeV/nucleon|iil, for it~ and K+, and at 10.8 GeV/nucleonU, for it~ . 
Lines show Gaussian fits to the sources. 





Ro [fm] 


A 


A(35fm) 


K+ (11.4 GeV/A) 


2.76 


0.702 


0.86±0.56 


TT- (11.4 GeV/A) 


6.42 


0.384 


0.44±0.17 


TT- (10.8 GeV/A) 


6.43 


0.486 


0.59±0.22 



Table 3: Parameters of Gaussian fits to the sources and integrals over imaged regions for 
the central Au + Au reactions. 



of emission pointsJor kaons and pions were predicted long ago within RQMD 
by Sullivan et al. E3. In the model, they were able to separate the different 
contributions to the source functions. 

The effects of long-lived resonances, mentioned above, are apparent in the 
sources extracted from the data. Thus, Table || gives the source integrals over 
the imaged regions together with parameters of the Gaussian fits to the sources. 

The errors are quite small for the fitted values. We find A^- < \k+ < 1 
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and A < A. 



7 Conclusions 

We have demonstrated that a model-independent imaging of reactions is pos- 
sible. Specifically, we have carried out one-dimensional imaging of pion, kaon, 
proton, and IMF sources. The three-dimensional imaging of pion sources is 
in progress. Our method of optimized discretization allows us to investigate 
the sources on a logarithmic scale up to large distances. The sources generally 
contain information on frcczc-out phase-space density, entropy, spatial den- 
sity, lifetime and size of the freeze-out region, as well as on resonance decays. 
The imaging gives us access to the spatial structure required to extract that 
information. 
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